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Abstract 

In this paper, a comprehensive spatio-spectral and temporal analysis for Chromobacterium violaceum colonies is 
reported. A hyperspectral imaging (HSI) system is used to recover the spectral signatures of pigment production in a 
non-homogeneous media with high spectral resolution and high sensitivity in vivo, without destructing the sample. 
This non-contact sensing technique opens avenues to study the temporal growing of a specific section in the 
bacterial colony. Further, from a 580 [nm] and 764 [nm] spatio-spectral time series, a wild-type and mutant 
Chromobacterium violaceum strains are characterized. Such study provides quantitative information about kinetic 
parameters of pigment production and bacterial growing. 
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Introduction 

For decades, natural pigments have been extensively used 
in various fields of everyday life such as food production, 
textile industries, paper production, agricultural practices 
and researches, water science and technology (Arad and 
Yaron 1992; Sirimanne et al 2006). Natural pigments not 
only have the capacity to increase the marketability of 
products, they also display advantageous biological activ- 
ities as antioxidants and anticancer agents (Stahl and Sies 
2005; Kong et al 2010). Several intensely colored com- 
pounds have been isolated from certain microorganisms 
and from various environmental sources (Ahmad et al. 
2012). An example is the violacein pigment, a purple- 
colored dye produced by one of the strains of Chromobac- 
terium violaceum, which is an indole derivative (Duran 
and Menck 2001) with antitumoral, antibacterial, antiul- 
cerogenic, antileishmanial, and antiviral activities. (Leon 
et al. 2001; Duran et al. 2007). 

The detection of bacterial components like antibi- 
otics, enzymes and secondary metabolites like pigments, 



"Correspondence: mgallardo@cefop.udec.cl 

1 Center for Optics and Photonics, Universidad de Concepcion, Concepcion, 
Chile 

Full list of author information is available at the end of the article 



is important in biotechnological industries (Steele and 
Stowers 1991). Traditional methods are effective but can 
take hours or days to realize them. Besides, they are 
destructive as well as time consuming. On the other hand, 
molecular detection based techniques for bacteria iden- 
tification are rapid, specific and sensitive. However, most 
are technically complicated and costly, and require well- 
trained specialists (Steele and Stowers 1991; Lazcka et al. 
2007; Velusamy et al. 2010). Now, molecular spectroscopy 
provides less invasive or even non-destructive methods 
that have been proven well-suited to investigate the chem- 
ical composition of biological systems (Cen and He 2007). 

Hyperspectral imaging (HSI) is an emerging technique 
that assimilates spectroscopy and imaging to provide both 
spectral and spatial information of a biological sample 
(Plaza et al. 2009). Further, the potential that HSIs have 
for in vivo optical diagnostics have being exploited due to 
the non-invasive feature and the massive spatio-spectral 
information collected by the system ( Vo-Dinh 2004; Pisani 
et al. 2013). In this sense, HSI systems record, for each spa- 
tial location being imaged, a set of hundreds of high spec- 
tral resolution images that jointly conform the spectrum 
of the biological sample. Such collections of images is 
known as hypercube (Borengasser et al. 2008), from which 
the pertinent qualitative and quantitative information is 
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Figure 1 Spot plating and simplified reflection based 
measurement scheme. Chromobacterium violaceum wild-type strain 
ATCC 31 532 (purple spots) and C. violoceum CV026 (white spots) are 
inoculated on a petri dish in LB agar. The plate is at uniform room 
temperature along the experiment. 



obtained (Sankaran et al 2010; Siripatrawan et al 2011). 
Different optical techniques can be applied to generate 
the hypercube from biological samples, some of them 
are based on tunable optical filters (Gat 2000), imaging 
spectrometers (Jun et al 2009), and coded hyperspec- 
tral imaging (Studer et al 2012), among other. In par- 
ticular, push-broom hyperspectral cameras (PBHCs) are 
imaging devices that carried out the spectral decompo- 
sition by means of an optical process (Kim et al. 2010). 
Further, PBHCs employ a scanning procedure to record 
spectral information of one spatial line at a time, with 
the ability to scan multiple batches of samples simul- 
taneously by moving across the process line. Also, the 
spatial information is important for monitoring the sam- 
ple as it can be used to extract its chemical mapping 



at different spatial points (Zavattini et al. 2003; Singh 
et al. 2010). In this regard, we can use HSI system to 
measure the spectral signature of chemical components, 
like bacterial pigments, inside the cells without altering 
the sample and its measurements even in a heteroge- 
neous medium where spectrophotometric signal is not 
resolutive. 

In this work, a spectral characterization of Chromobac- 
terium violaceum pigment through a HSI is performed 
from the visible spectrum to part of the near infrared spec- 
trum (400 nm to 1000 [nm] with 1,04 [nm] of spectral 
resolution). It is well known that a pigment changes the 
color of the reflected light as a result of selective color 
absorption. In this sense, the pigment is quantified associ- 
ating its signal to the maximal optical absorption centered 
at 580 [nm] (Gerhardt, 1970; Koch, 1994). With HSI it 
is now possible to resolve signals and quantify the bac- 
terial pigment with a remote, non-invasive and real-time 
procedure comparatively to a spectrophotometric mea- 
surement. Finally, a series of spectral signatures masking 
relevant information that can serve as an indicator of a 
specific bio-processes is obtained, i.e. pigment production 
and bacterial growing rate. In this case, despite using a 
small number of spectral bands, these are selected as a 
function of the amount of spectro-temporal information 
associated with the pigmentation process. For a differ- 
ent process, the evaluation of the spectral information 
can determine a larger number of spectral bands required 
for characterization. Hence, the amount of useful spec- 
tral information is directly related to the process being 
measured. 

Materials and methods 

Bacterial cultures 

The C. violaceum wild- type strain ATCC 31532 and C. 
violaceum CV026 (mini-Tn5 mutant of ATCC 31532) 
(McClean et al. 1997) are cultured in Luria Bertani (LB) 
broth (10% peptone, 5% yeast extract, 5% NaCl, per 1L 
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Figure 2 Laboratory prototype to retrieve the spatio-spectral image. The aluminum arc holds the illumination and camera system. The target 
scene lies on a movable base where the arrow indicates scan direction. Imaging is acquired during the scan process and data analysis is 
post-processed. 
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Figure 3 Selected region of interest in the reasembled spectral image at 580 [nm] at a) t = 0 hrs and b) f = 37 hrs. 



distilled water) (Sambrook et al. 1989). Solid bacterial 
medium is made by the addition of agar at a final concen- 
tration of 15 g/L. 



measurement an aliquot of liquid culture is taken and the 
absorbance is determined every 50 [nm] (between 350 and 
700 [nm]). 



Spectrophotometric measurement 

The bacterial population and the pigment production 
is estimated by measuring absorbance of the culture 
using a spectro-photometer (UV-1700 UV-VIS Spectro- 
photometer) at different wavelengths (350-700 [nm]). 
Measurements are taken each one hour during lag phase 
and every 30 minutes in the exponential phase. For each 



Spot plating and incubation 

The initial liquid bacterial culture is grown overnight 
at 30°C in a thermo-regulated bath. Agar LB plates are 
inoculated with 20 fiL spots at pre-determined and well- 
spaced locations on the agar surface. We separated the 
wild-type strain from the mutant strain to avoid con- 
tamination. Figure 1 shows a scheme of the imaged spot 
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Figure 4 Block diagram of the reflectance recovery process. The sequence start with the scanning of the target scene. Data are then 
post-processed, here the calibration model compensates for the system spectral response which gives the spatio-spectral data cube. Then, the 
spatio-spectral reading identifies the specific wavelength of maximum absorbance that allows to characterize the pigment intensity change. Finally, 
the pigment signal is plotted as a function of time. Here, each point is a scanning sequence at that time. The scanning frequency is of the order of 
hours for the lag and stationary phase and half an hour for the exponential phase. 



Gallardo etal. AMB Express 2014, 4:4 
http://www.amb-express.eom/content/4/1/4 



Page 4 of 9 



locations. The plate is incubated at room temperature in 
the camera stage for 24 hours. Three measurements are 
made, each one with four replicates of the wild-type strain 
and two replicates of the mutant strain. 

Experimental Setup 

In order to acquire image sequences with high spectral 
resolution, The set-up illustrated in Figure 2 is con- 
structed, where the key component is the PBHC. Such a 
camera can be described, in brief, as an optoelectronic 
imaging system composed of four sub-systems: the optics, 
the spectrograph, the sensor, and the readout electron- 
ics (Borengasser et al. 2008). Unlike broadband images, in 
PBHC the spatial information is mapped onto one axis of 
the focal plane array (FPA), while the spectral information 
is mapped onto the second axis. Therefore, a biological 
sample (now target scene) must be scanned one line at 
a time by displacing the camera in an orthogonal trajec- 
tory with respect to the axis used to encode the spectral 
information. 

The HSI system employed in this work is the Photon- 
focus Hurricane 40 V10E PBHC. This camera is based on 
a CMOS FPA of 1024 x 1024 photo-detectors. Further, 
the input radiance is decomposed and converted into sev- 
eral spatio-spectral images through an optical process that 
allows us to obtain spectral bands between 400-1000 [nm], 
with a spectral resolution of 1.04 [nm]. Also, due to the 
nature of the acquisition process of PBHC, a mobile plat- 
form has been constructed to simulate and synchronize 
the scanning procedure. 

Furthermore, to compensate for any noise source 
degrading the hyperspectral data, it is required to perform 
calibration measurements. To this end, we have employed 
a quartz tungsten halogen (QTH) calibration lamp, which 
guarantees an uniform and continuous spectral illumi- 
nation between 200-2500 [nm]. Also, a Spectralon SRT- 
99-120 is used as a diffuse reflectance target that ensures 
a reflectance of 99% between 250-2000 [nm]. 

As shown in Figure 2, the diffuse target lies on a 
petri dish holder to increase the reflectivity of the semi- 
transparent agar and colonies. Further, during the scan- 
ning process, each line of the target scene acquired 
receives the same angle of illumination, then we can 
ensure that all spatial points were illuminated under the 
same conditions. The sample measurement scheme is dis- 
played in Figure 1, where the input radiance is generated 
by the light reflected from the agar. Finally, the dimensions 
of the hypercubes produced after scanning the whole tar- 
get scene are 1024 spatial pixels, 574 spectral bands, and 
1000 temporal lines. 

Spectral reflectance recovery 

In order to retrieve the spectral reflectance of biological 
samples, several stages have to be performed. First, the 



sequence starts with the scene scanning due to the nature 
of the of the PBHC acquisition process. Second, the input 
radiance of a spatial line is decomposed and sampled into 
P spatial pixels and Q spectral bands, taking S line sam- 
ples at different times. Hence, the entire target scene is 
mapped into a hypercube of dimensions Px QxS. Third, 
the hyperspectral data must be calibrated to compensate 
for any degrading effect produced during the signal trans- 
duction. To do so, we have mathematically represented 
the hyperspectral response by the following first-order 
model: 

Y{U k h k) = r(kj)a(i,j)R(i, X Jf k)+b(i,j) + V(i,j, k), (1) 

where the suffixes i, kj, and k denote, respectively, the 
spatial location at the scene, the spectral band, and the 
temporal sample, where i = 1, . . . ,P, j = 1, . . . , Q r , and 
k = 1, . . . , S. The variable R(i, Xj, k) represents the input 
reflectance collected at the kjth. spectral band by the ijth 
photo-detector at the kth sample time. The terms a(i,j) 
and b(i,j) represent, respectively, all the multiplicative and 
additive noise sources corrupting the output of a PBHC. 
The additive term V(i,j,k) is known as the temporal 



a 



3 

d 

■e o. 

o 

in 

3° 



□ 


350 nm 


o 


400 nm 


▲ 


450 nm 


V 


475 nm 


o 


500 nm 


A 


525 nm 


T 


550 nm 


□ 


600 nm 


▲ 


700 nm 



■ A A 

i I ! 



4 

Time(h) 



3 

d 

■e o. 

o 

m 

3° 



□ 


350 nm 


o 


400 nm 


▲ 


450 nm 


V 


475 nm 


o 


500 nm 


A 


525 nm 


T 


550 nm 


□ 


600 nm 


A 


700 nm 



6 

A 

■ i ¥ 



1 



1 if 



2 4 6 8 

Time(h) 

Figure 5 Spectrophotometric measurement of both C. violaceum 
wild-type strain (a) and C. violaceum CV026 strain (b) in liquid 
culture. Data shows the absorbance intensity signal versus time in 
the visible spectrum. 
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Figure 6 Spectrophotometric measurement of both C. violaceum 
wild-type strain and C. violaceum CV026 strain in liquid culture. 

Data shows the absorbance intensity signal versus wavelength. Data 
is taken during the stationary phase of growing, seven hours after 
inoculation. 



white noise associated with the readout electronic for the 
ijxh photo-detector. However, in most of the cameras the 
white noise V(i,j,k) is negligible compared to the term 
b(Uj), so it can be disregarded. Finally, the parameter r(Xj) 
represents the system spectral response, i.e. the slit, the 
diffraction grid response and the spectral responsivity of 
the detector material. 

In general, the additive noise term b(i,j) is treated as 
a baseline value produced by different noise sources, i.e., 
dark current, non-uniformity and amplifiers noise. Such 
baseline can be estimated through the measurement of 
a dark reference. This is, considering that the spectral 
decomposition is carried out by an optical process, block- 
ing the camera input prohibits the decomposition process, 
hence, R(i, kj, k) = 0. Further, both multiplicative terms, 
r(Xj) and a(i,j) can be estimated by means of a white ref- 
erence, which is obtained by illuminating a diffuse target 



with the QTH calibration lamp. Consequently, the fol- 
lowing formula is used to compensate and retrieve the 
spectral reflectance: 



R{U A./, k) 



Y(i,Xj,k)-b(i,j) 
Y r (i,kj,k)-b(i,j) 



(2) 



where Y r (i, Xj, k) is the white reference value and b{i,j) is 
the estimated baseline obtained when measuring the dark 
reference. Recalling Eq. 1, by subtracting b (i,j) in Eq. 2, the 
baseline value is compensated. Further, the normalization 
between the measured sample and the white reference 
compensates for the system spectral response and the 
term a(i,j). With the spectral reflectance at hand, the 
absorbance values A(i,Xj,k) are retrieved and the maxi- 
mal and minimal pigment absorption bands are selected, 
that in our case are 580 [nm] and 764 [nm], respectively. 
Therefore, to analyze the spectral behavior of the samples 
taken at different times, first it is required to reassem- 
ble the spatial information from the hypercube. Then, a 
region of interest region of interest (RI) is located on each 
sample, covering a spatial region equal to 21 by 21 pixels, 
as shown in Figure 3. 

Each absorbance value enclosed by the RI is averaged to 
minimize any small variation produced by sensor noises 
or by simple irregularities. Thus, the sample averaged 
absorbance is calculated as follows: 

/o+10 /ro+10 

A(i 0 ,Xj,ko) = A (n,Xj,m) (3) 

n—io-10 m—k{)-\0 

where (io,ko) is the RI center. It should be noted that 
since each pixel represents an absorbance value, select- 
ing random measurement points is also a valid alternative. 
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Figure 7 Imaging at 580 nm of both C. violaceum wild-type strain (a) and C. violaceum CV026 strain (b). Hyperspectral measurements in 
(c) show the absorbance intensity signal versus time for both strains. 
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Figure 8 Imaging at 764 [nm] of both C. violaceum wild-type strain (a) and C. violaceum CV026 strain (b). Hyperspectral measurements in 
(c) show the absorbance intensity signal versus time for both strains. 



Finally, to summarize the methodology, a block diagram 
explaining each stage of the reflectance recovery process 
is presented in Figure 4. 

Results 

One of the advantages of using HSI systems is to recover 
the spectral signatures of pigment production in turbid 
media with high spectral resolution and sensitivity, in con- 
trast to optical absorption based instruments. The tradi- 
tional spectrophotometer procedure disposes the sample 
inside small quartz containers. Here, colored laser beams 
creates the observed absorption spectrum. When the 
sample is a complex mixture of bacterial cell components 
and pigment, the observed spectrum is indistinguishable 
for a clear identification, e.g., between the C. violaceum 
wild-type strain and C. violaceum CV026 strain as shown 
in Figure 5. 

Although both samples are differentiable trough sim- 
ple eye-inspection, the absorbance intensity signal shown 
in Figure 6 is poor for quantification. This constraint is 
due to the screening by the turbid media itself which is a 
result of optical absorption in the visible spectrum by the 
bacterial cell components (membranes, organelles, DNA, 
proteins, etc.) (Gerhardt 1994; Koch 1970). 

With the aim to resolve the spectral differences between 
the strains associated with the pigment production, we 
use a PBHCs that allows to obtain spatio-spectral infor- 
mation without disturbing the sample. The experiment 
is developed by plating both strains in LB-agar at equal 
initial cell concentration, and following the methodol- 
ogy explained in the previous section. In particular, two 
spectral bands are selected and the subsequent temporal 
samples have information regarding the biological pro- 
cess. To obtain the signal, a band centered at 580 [nm] 
which corresponds to the wavelength of maximal absorp- 
tion, is selected. In Figure 7c is shown the absorbance 



intensity signal versus time for both strains. Now 
the pigment absorbance signal is successfully retrieved 
and is clearly differentiable between the strains. The 
absorbance pigment intensity signal of the C. violaceum 
wild-type strain exhibits a bacterial growing-like curve 
behavior. 

In order to obtain quantitative information regarding 
the pigment production curves of Figure 7, some consid- 
erations need to be taken into account. The absorbance 
at 580 [nm] has simultaneous information of pigment and 
bacterial cell components. Thus, a band centered at 764 
[nm] which corresponds to the minimal pigment absorp- 
tion, is now selected. The absorbance intensity signal 
versus time for both strains at 764 [nm] can be seen in 
Figure 8. 

The information regarding pigment and cell compo- 
nents can be summarized in Figure 9. Data is fitted by 



wild-type 580 [nm 
mutant 580 [nm] 
wild-type 764 [nm] 
mutant 764 [nm] 




15 20 25 
Time (h) 

Figure 9 Fitting of hyperspectral measurements of both C. 
violaceum wild-type strain and C. violaceum CV026 strain at 
580 [nm] and 764 [nm]. 
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Table 1 Kinetic parameter of pigment production and bacterial growth 



Hyperspectral imaging measurement 3 



C. violaceum 
wild-type strain 



C. violaceum 
CV026 strain 





g(h) c 


kC/?" 1 ) 0 


g(h) c 


kc/?- 1 ) 0 


580 [nm] (pigment production + bacterial growth) 


5.28 ± 1.10 


0.1 37 ±0.029 


29.77 ± 6.85 


0.025 ± 0.006 


764 [nm] (bacterial growth) 


7.57 ± 1.14 


0.094 ±0.014 


26.93 ±7.31 


0.028 ± 0.008 


580 [nm]/764 [nm] (pigment production) 


15.9 


0.043 


d 


d 


OD measurement 


2.93 ± 0.65 


0.249 ± 0.055 


8.4 ± 1 .23 


0.08 ± 0.01 1 



a measurements are made on solid medium with no agitation at RT. 
Measurements are made on liquid medium with agitation at RT. 

c The coefficients are estimated considering a confidence level of 95% that the real value is in that range. 
d mutant strain does not show spectral response. 



using a non-linear least squares method on the exponen- 
tial phase. Kinetic parameter of this fitting are shown in 
Table 1. 

The grown parameters g and k for both strains at 764 
[nm] are calculated. At this wavelength only information 
of bacterial growing is obtained. Although the wild-type 
strain grows faster than the mutant strain both have 
comparable velocities k. This is in agreement with the 
values found in literature and those obtained with the 
spectrophotometric measurements (see Table 1). 

To quantify the violacein production, a division of func- 
tions is obtained. The foregoing is achieved between 
the curve regarding information of the bacterial pig- 
ment and bacterial growth (Figure 9) -squares and the 
curve containing information only of the bacterial growth 
(Figure 9)-triangles (Figure 10). Comparing the veloci- 
ties of mutant strains at 580 [nm] and 764 [nm], noted 
that both values are equal, as expected, due to the lack 
of pigment production. C. violaceum CV026 strain are 
does not show spectral response in the spectral range 
analyzed. 



Discussion 

The violacein production can be quantified using a 
spectrophotometer, after an extraction and purification 
process (Rettori and Duran 1998; Lu et al. 2009). This 
procedure fails to identify pigment remotely and not 
allows to follow the process in time. To follow the pig- 
ment production without destroying the bacterial cells we 
use a Hyperspectral Imaging System. This system provide 
spectral and spatial information of the Chromobacterium 
violaceum enabling us to obtain kinetic parameters of 
bacterial growth and pigment production. 

The violacein production can be quantified at 580 [nm] 
by subtracting the contribution of the bacterial cell com- 
ponents at 764 [nm] which resulted in k p i gment = 0.043. 
Therefore, the results are in good agreement with the 
values calculated using optical density (OD) at 600 [nm] 
(see Table 1) and with those found in the literature 
(Corpe 1964). 

In this work, a real time, non invasive and simple process 
to quantify bacterial pigment production trough HSI sys- 
tem is reported. The PBHC camera produces a hypercube 



5 




wild type 580/764 

mutant 580/764 

0 I i 1 i i 1 1 

0 2 4 6 8 10 12 

Time (h) 

Figure 10 Mathematical division of the fitting curves generated at 580 [nm] and 764 [nm] of both C. violaceum wild-type strain and C. 
violaceum CV026 strain. 
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with spatial and spectral information in function of time 
without destroying the sample. Hence we can analyze 
multiple samples with high wavelength resolution in a sin- 
gle experiment. Based on the process, the technique can 
be extensible to multiple colored samples at the same time. 
This optical technique could be useful for the analysis of 
complex samples such as bacterial biofilms, environmen- 
tal samples or non-culture bacteria (Cortizo and de Mele 
2000). This can be applied to quantify another biologi- 
cal pigments obtained from natural source without alter 
the environmental conditions. For example, studies can 
be done to quantify photosynthetic pigments, like chloro- 
phylls or carotenoids, presents in some plants or bacteria. 
Also this technique can be used in the pigment determi- 
nation in fruits or other exportation products. The advan- 
tage of this technique is the decomposition of the compo- 
nents of the biological system in a set of wavelengths each 
one associated to a specific spatio-temporal character. 
This allows to get new insights of a particular process in 
a complex media. Also, since we have detected the spe- 
cific wavelengths that allows to characterize the bacterial 
pigment production and growing process, the next step is 
developing a filter-based imaging system with the aim to 
increases the temporal resolution of the measurements. 
That is, designing a monochromatic camera that includes 
a filter tuned at the wavelength of interest, discarding the 
scanning stage and acquiring the spatio-spectral image 
that exhibits greater variation in time (Nishino et al. 
2013). Further, such devices do not require a scanning 
stage, which could increase the temporal resolution of our 
measurements. 
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